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A Search for Neutron Star Precession and Interstellar Magnetic 
Field Variations via Multiepoch Pulsar Polarimetry 
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ABSTRACT 

In order to study precession and interstellar magnetic field variations, we 
measured the polarized position angle of 81 pulsars at several-month intervals 
for four years. We show that the uncertainties in a single-epoch measurement 
of position angle is usually dominated by random pulse-to-pulse jitter of the 
polarized subpulses. Even with these uncertainties, we find that the position 
angle variations in 19 pulsars are significantly better fitted (at the 3a level) by 
a sinusoid than by a constant. Such variations could be caused by precession, 
which would then indicate periods of ~ (200 — 1300) d and amplitudes of ~ 
(1 — 12) degrees. We narrow this collection to four pulsars that show the most 
convincing evidence of sinusoidal variation in position angle. Also, in a handful 
of pulsars, single discrepant position angle measurements are observed which 
may result from the line of sight passing across a discrete ionized, magnetized 
structure. We calculate the standard deviation of position angle measurements 
from the mean for each pulsar, and relate these to limits on precession and 
interstellar magnetic field variations. 

Subject headings: ISM: magnetic fields — polarization — pulsars 

1. Introduction 

Many pulsars are among the most highly polarized astrophysical sources of radiation. 
One polarized parameter, the position angle of linear polarization tp (also called the PPA - 
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polarized position angle) is a particularly useful quantity: it carries information about the 
i ntrinsic geometry of the emitt ing; region at the pulsar [e.g., via the rotating vector model 
( IRadhakrishnan fc Cookdll969l )]; it also probes the magnetized plasm a along the line of sight 
through Faraday rotation measurements [e.g., IWeisberg et al.l (J2004J)]. The purpose of this 
paper is to study the behavior of the polarized position angle on timescales of 10 2 ~ 3 d, in 
order to investigate time-variable behaviors of pulsars' magnetospheric geometry and the 
interstellar medium. For example, precession of the neutron star or other long-timescale 
oscillatory behavior affecting the emission beam will change the PPA, as will variations in 
the magnetized interstellar plasma along the line of sight. 



2. Observations and Analyses 

Ninety-eight pulsars were observed at 21 cm from Arecibo Observatory every few months 
for several years from 1989 to 1 9 93. T he resulting grand average polarized pulse profiles 



were analyzed in IWeisberg et al.l ( 119991 ). where observing details can also be found. In the 
present work, we found that we had adequate data on eighty-one of the pulsars to search 
for temporal variations of polarized position angle. For all of them, the polarized position 
angles were originally measured relative to an unknown origin. Consequently, in the present 
work, in order to search for temporal changes in position angle, we needed to reference our 
data to a position-angle standard. To do so, we selected two calibrator pulsars, B0540+23 
(J0543+2329) and B1929+10 (J1932+1059). The calibrators were chosen on the basis of 
the relative constancy of their PPA's during full-sky Arecibo observing tracks. We then 
compared the measured polarized position angle of a target pulsar on a given observing 
session to the PPA of one of the two calibrator pulsars. While the absolute position angle is 
still unknown, this referencing procedure enables us to be sensitive to position angle changes 
over time, assuming only that the position angle of the calibrators remains fixed over time. 
The procedures are described in the following sections; mathematical details are given in 
Appendix A. 

In the rare cases where a calibrator pulsar was not observed during a daily session, we 
created an artificial "meshed" calibrator pulsar by determining its expected position angle 
on that session from observations of all other pulsars during that session. 
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2.1. A Pulsar's Polarized Position Angle Curve 



Pulsars tend to be very highly linearly polarized objects. The observed linear polariza- 
tion or iginates at or near the star, of ten following the magnetic field lines at the point of 
origin (IRadhakrishnan fc Cookelll969f ). and is then modified by various propagation effects. 
We define the generic polarized position angle "curve " for a pulsar at epoch t, a function of 
pulsar longitude at fixed t, as the set of polarized position angles i/j(4>,t) for < < 2n. 
(See the bottom panels of Fig. [2] for examples of PPA curves.) 

The measured PPA curve at epoch t can be expressed as the sum of five terms - an 
ensemble average (denoted below by angular brackets) for the indicated pulsar longitude 
and epoch t which reflects the pulsar's intrinsic orientation, an interstellar Faraday rotation 
term, an instrumental origin, and a term arising from errors due to radiometer noise and one 
from intrins ic pulse-to-pul se fluctuations in amplitude and phase of individual pulses (called 
"jitter;" see ICordesI (119931 ) and further discussion below): 



■0(0, t) = (0(0, t)) + 4>Faraday(t) + l/j instr (t) + 5^ noise ((f), t) + 5ljj jitter(4> ', t). (1) 

The first two terms on the right-hand side of the above equation hold the physics of interest in 
this experiment - namely the temporal behavior of the pulsar's orientation or magnetosphere, 
and of the interstellar magnetic field through which its signals propagate, respectively. The 
third term defines the instrumentally measured PPA corresponding to an absolute PPA on 
the sky of zero. The latter two terms, though having zero mean, represent phenomena that 
add variance to the measurements, thereby limiting our precision in specifying the first two. 
We now discuss each of these five terms in order. 



2.1.1. Variations in PPA Due to Changes in Orientation of the Pulsar or its 

Magnetosphere 

(a) Precession: Precession would lead to a periodic wobble of the pulsar orientation. 
Most pulsar magnetospheric models indicate that the polarized emission is tied fundamen- 
tally to the geometry of the magnetic field in the emitting regions. Hence it is expected that 
a periodic wobble in pulsar orientation translates directly into a periodic variation of similar 
amplitude in polarized position angle. Therefore, a search for temporal variations in PPA is 
a sensitive probe for the presence of precession. 

The simplest model is one where the wobble would result from the free precession of an 
asymmetric isolated pulsar, although there also exist models where the precession is forced 
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by an external disk or companion (jQiao et al.ll2003f ). or by the star's own radiation (IMelatos 



20001 ) . We assume that the neutron star is approximately spherical, with a slight axially 



symmetric oblateness about a symmetry axis. Hence the three moments of inertia have the 
relationship I x = L\ > 1% — J 3 , where I\\ is the moment pa rallel to the symmetry axis. 
(For e xt ension of the p reces sion model to a triaxial bo dy, see lAkgiin. Link, fc Wasserman 
( 120061 ) ) . IShahaml (119861 ) and iJones fc Anderssonl (120011 ) show that the motion of a point on 
the star can be decomposed into a quick rotation about the total angular momentum vector 
plus a slow precession about the body symmetry axis. For a sketch of the process and a 
link to an animation, see Fig. [TJ for readers wishing to modify the geometrical parameters 
of the a nimation, the Mathemati ca source code is also available online at the Mathsource 
library ( jEverett fc Weisbergl 120071 ). The angle 9 between the angular momentum vector and 
symmetry axis leads to a "wobble" at both the rotation and precession frequencies, with the 
latter leading directly to a change in PPA of similar amplitude^ on precession timescales. 
Very slow precession, with a period of many years, would produce just a linear trend in PPA 
with epoch. 

(b)Magnetospheric Adjustments: As we will describe below, there are theoretical 
difficulties (though not insoluble ones) with the picture of precession-induced pulse shape 
variations, c entered primarily on the long observed timescales and relatively large inferred 
amplitudes. iRuderman fc Gill ( 120061 ) have proposed an alternate model whereby the pulse 
shape and timing variations are caused by slowly rotating magnetospheric currents rather 
than by precession of the whole star. While this is an interesting alternative interpretation 
of the observations, more work needs to be done to flesh out the model and to verify that it 
can explain the observed timescales and amplitudes. 



2.1.2. Variations in PPA Due to Changes in Interstellar Faraday Rotation 

In this section, we focus on studying variations in the magnetized ISM plasma via a 
search for temporal variations in PPA. The second term of Eq. [T] quantifies Faraday rotation 
of the PPA in the interstellar medium, ip Faraday {t), at some epoch, t. The rotation measure, 
RM(t), is defined so as to specify the angle of PPA Faraday rotation at wavelength A: 

^Faraday(t) = \ 2 RM(t), (2) 

with angle in radians, A in m, and RM(t) in rad mT 2 . 



1 The amplitude of the PPA variation, Aip, is magnified by a factor of ~ [cos a] 1 ; .i.e, ~ 9 x [cos a] ' 
where a is the (generally unknown) latitude of the magnetic axis with respect to the symmetry equator. 
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A measured change in PPA, AipFaraday, then implies a change in Faraday rotation mea- 
sure, ARM: 

Alp Faraday 



ARM 



X 2 



(3) 



In order to relate ARM to the ISM, we express RM as a function of fundamental 
properties along the pulsar-Earth path: 



RM = 0.81 



D 



Bu 



,s n e \s 



ds, 



(4) 



where B\\ is the magnetic field parallel to the line of sight in /iG, n e is the electron den- 
sity in electrons cm" 3 , and ds is a differential path element whose integrated distance is D 
( iLyne fc Graham-Smithl 120061 1 . 



Each of the variables in Eq. H] may be a function of time. While this time- variability 
could result directly from temporal changes in the ISM, it would more likely be induced by the 
pulsar's motion rapidly sweeping the line of sight across a presumed spatially inhomogeneous 
ISM. 



Since measurements show that temporal variations of dispersion measure DM 



are generally quite small (IPhillips fc Wolszczan 



1991 



Weisberg et al.l 12004 lAhuja. Gupta, fc Kembhavi 



Backer et al. 



2005 



You et al. 



1993; Hobbs et al.l 2004b 



F_ rtp ds 



20071 ) . we further expect 



that any observed variations in RM would originate principally from the pulsar carrying the 
line of sight across an inhomogeneity in Si. Therefore, we model the RM change as resulting 
from the line of sight passing across a discrete region having a variation in the parallel field, 
ABu, occupying a fraction / of the pulsar-Earth path, with no associated n e variation. 



In this case, ARM, the RM change due to the discrete region is 

rfD 

ARM = 0.81 / AS|| n e ds « 0.81 AB\\ f DM. 
Jo 



(5) 



Finally, we can relate an observed change in PPA, Aip Fara day, to the properties of the 
discrete region by combining Eqs. |3] and [5j 



/ AS., 



ARM Aip Faraday 



0.81 DM 0.81 A 2 DM' 
or, expressed clS Si fractional variation: 

/ A£|| _ ARM = A^ Farad g y 
(£,[) ~ RM \ 2 RM ' 

where (Bu) is the line of sight average. 



(6) 



(7) 



-6 - 



2.1.3. The Instrumental PPA and its Variance 

The third term in Eq. [I] represents a conversion of absolute PPA on the sky, to PPA 
measured in the instrumental frame. The instrumental frame involves geometrical terms 
associated with telescope orientation as well as instrumental phase delays during a given 
session. As noted above in §2J we did not directly determine the instrumental PPA but 
rather referenced our measured PPA to the PPA of a calibrator pulsar, which has its own 
variance (see §2.2p . 

2.1.4- The Variance in PPA due to Radiometer Noise 

The noise term in Eq. [T] results from the finite sampling statistics of a noisy signal and 
is familiar to all radioastronomers. Its contribution to the variance, noise , was determined 
empirically by measuring the variance off-pulse. 

2.1.5. The Variance in PPA due to Pulse Jitter 

The jitter term in Eq. [T]is unique to pulsars and requires further explanation, especially 
because it often dominates the error budget. We characterize as "jitter" all PPA changes 
on pulse-to-pulse timescales that are counterparts to the arrival-time jitter of individual 
subpulses. These random PPA variations occur due to various propagation effects in the 
pulsar magnetosphere analogous to Faraday rotation. These fluctuations might be correlated, 
for example, for the duration of a single subpulse in a given pulse period, but would be 
uncorrelated between subpulses and across pulse periods. Random variations in altitude 
of emission would produce a PPA curve that is shifted earlier or later for higher or lower 
altitudes, respectively. This kind of variation would also be highly correlated across a single 
subpulse. The variance of ip due to jitter, j itter , is calculated as follows. 

By definition, uncertainties <T/ j jitter (0) and 0£, j itter (0) in polarized parameters I ((f)) (to- 
tal power at longitude bin 0) and L (linearly polarized power at longitude bin 0) due to 
random pulse-to-pulse jitter are characterized by 

07, jitter (0) = mi(4>)I((j)), (8) 

and 

<7L, jitter (0) = m L ((j))L((j)), (9) 

where mx is the "modulation index" of polarized parameter X, a measure of its pulse-to- 
pulse modulation. 
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In our experiment, each polarized profile was integrated online for two minutes, thereby 
attenuating these pulse-to-pulse modulations. Therefore, the jitter-induced fluctuations in 
the measurement of these quantities after N pulses, will be 



07, jitter — jitter ((/>)/ VN = mi {(f>)I(<f>)/VN, (10) 

and 

jitter — jitter (4,)/VN = m L {<t>)L(<t>)/VN. (11) 

At high S/N and under the reasonable assumption that ctq ~ au, o L propagates to an 
uncertainty in position angle 

^(0) « a L {<f>)/{2L{<f>)) (12) 



(jNaghizadeh-Khouei fc Clarkd 119931 ). Therefore, the jitter-induced uncertainty in the mea- 



surement of PPA in an integration of iV pulses is 

C"V Jitter {<t>,N)&m L {<t>)/{2VN) (13) 



As we shall see below in §2.21 we must ultimately determine critter, the jitter-induced 
fluctuation in the PPA integrated not only over N pulses but also over the full range of 
longitudes. With the subpulse as the basic unit of correlated emission, the jitter-induced 
fluctuations in the measurement of PPA after N pulses, each containing n subpulses, will be 
characterized by 

Ctyjitter = 0"V jitter (N, n) = (Tetter (0, N) / y/n PS m L / (2\fNn). (14) 



Since we could not directly measure the quantities tul and n with our two-minute 
integrati ons, we surveyed t he literature for typical values. Many observ ations suggest that 
mi ~ 1 ( Bartel et al. 1980 : Weisberg et al. 1986 : Weltevrede et al. 20061) and the f ew extant 
polarized modulation index observations (e.g., IJenet. Anderson, fc Prince! (120011 )) support 
the reasonable assumption that mi ~ 1 as well. Single-pulse observations indicate that n is 
also approximately 1. Finally, we found that constant-?/' fits to the full set of pulsars had a 
mean xl ~ 1 when mi/s/n ~ 1.7 in Eq. [TU so we adopted this value; i.e., 



^jitter ~ 0.85/v^. (15) 

We empirically tested the relative contributions of radiometer noise and pulse jitter to 
the PPA uncertainty in the following tests on one of the two calibrator pulsars, B0540+23. 
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We observed it for several ~two-hour stretches with a series of standard two-minute integra- 
tions identical in form to all others in this experiment, and analyzed the resulting data as 
follows. First, a linear trend in PPA, presumably due to variable ionospheric Faraday rota- 
tion, was fitted and subtracted out [§. For PSR B0540+23, the uncertainty in PPA resulting 
from radiometer noise was Uip^noise ~ 0?2, while the observed standard deviation of two-minute 
integrations from the linear trend was a far larger ~ 1°. Pulse jitter considerations, however, 
lead via Eq. [15] to an estimated critter ~ 2°. 

While the jitter uncertainty estimated via Eq. [[^appears to be about twice the observed 
standard deviation, it is important to bear in mind that the simplifying assumptions used in 
the jitter calculation can easily lead to factor of ~ 2 errors in either direction for a particular 
pulsar, as they apparently have in this case. (See further discussion in §3.2p . Nevertheless, 
the important conclusion of this exercise is that the actual PPA uncertainty is far larger 
than that due to radiometer noise alone, and that pulse jitter appears to be the predominant 
source of the observed noise, with an amplitude in rough agreement with our model. Note 
also that the jitter model should give excellent results for a given pulsar if all factors in Eq. 
dHare directly measured, rather than using an average over our pulsar sample, as in Eq. [Tj3 



2.2. Estimation of a Pulsar's Polarized Position Angle Offset with Respect to 

a Calibrator Pulsar, and its Variance 

In order to study temporal variations in a pulsar's orientation or in the intervening 
medium, we investigate the PPA as a function of time. As noted in £J2] , the instrumental 
PPA origin, VWtr(^) of Eq. [TJ was not explicitly measured. Instead, we reference our target 
pulsar PPA's to "calibrator" pulsars B0540+23 or B1929+10, which can be considered PPA 
"beacons" of constant but unknown PPA. We do so in a three-step process: first we find the 
offset of the target p ulsar's PPA cury e at ep och t from a high S/N "template" PPA curve of 
the same pulsar from lWeisberg et al.l (119991 ) (see Fig. [2^ and Eq. IA1[) . Next we do the same 
procedure at the same epoch on the calibrator pulsar B0540 or B1929+10 (see Fig. [2b and 
Eq. IA2j) . Finally, to yield the desired offset of the target pulsar from the calibrator pulsar 



at epoch t, Aip(t), we difference the above two results (see Eq. IA3j) . 

The variance in a target pulsar-template PPA difference measurement, o^uy is the 
sum of the noise and jitter contributions from the target pulsar at epoch t (the template 



2 To first approximation, ionospheric Faraday rotation is nulled in our usual pulsar - calibrator pair 
measurements, but not in these all-sky measurements of a single pulsar. 
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contributes negligibly) : 

°%{t) = °%{t), noise + °%{t), jitter- (16) 

A similar equation holds for the variance in a calibrator pulsar-template PPA difference 
measurement, , t y However, we are ultimately interested in the variance of the position 
angle of a pulsar with respect to the calibrator pulsar, (cf. Eq. IA3j) . This variance, o^f*)' 
will then be the sum of the above-discussed variances in the pulsar-template and calibrator- 
template PPA difference measurements: 

°AW) = a m + °% cal (ty ( 17 ) 



3. Results 

We searched for evidence of temporal variations of polarized position angle in each of 
the 81 target pulsars, with an emphasis on a few likely signatures: (1) sinusoidal variations in 
PPA due to periodic phenomena at the pulsar such as precession, as described in §2.1.1( (2) 
discontinuous changes in PPA caused by interstellar propagation effects, such as the line of 
sight crossing a region with significantly different magnetic field, as discussed in §2.1.2} and 
(3) ramps in PPA caused by phenomena such as those listed above, having a characteristic 
time longer than our experiment's. For each pulsar, the dataset consists of the measured 
PPA offsets, Aip(t) (cf. Eq. IA3j) , and their estimated uncertainties, OAV>(t) ( c f Eq. [17)1 . at 
a set of epochs, t. In the following sections, we give an overview of these results. See the 
Discussion for exploration of the astrophysical implications. 



3.1. Sinusoidal Variations of Polarized Position Angle 

To test for possible periodic variations in each pulsar's PPA, we fit a sine function to 
Aip(t). This sine function requires four free parameters: the zero-point Aip about which 
the sinusoid oscillates; the amplitude of the sinusoid about Aip ; the period of the 
sinusoid; and a phase offset t : 

Ai/>(t) = A^o + H sin[27r(t - t )/P^}. (18) 

In fitting with these free parameters, we had two chief concerns: that our best-fit results for 
this function's parameters not be dependent on initial guesses for those parameters, and that 
we only accept the sinusoidal fit's results when that fit represents a statistically significant 
improvement over the fit of a constant, A^. 
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To ensure that the fit results were not dependent on initial guesses, we ran trial fits on 
a large number of initial values for both sinusoidal period and amplitude a^. The initial- 
parameter search space for the amplitude was set to lie in the interval [0°, 2 • max(A^(t) — 
Aip)\. Meanwhile, the maximum sinusoidal period was set to 5 • At, where At is the total 
timespan of the observations, in order to look for long-timescale variations. The minimum 
sinusoidal period tested was to 2 • At/A^ sess i on s where iV sess i ons is the number of observation 
sessions for each pulsar. This range of sinusoidal amplitude and period was then used to 
define a linearly-spaced grid of 501 x 501 initial guesses. A separate non-linear least-squares 
fit was run for each set of initial parameters, where all four parameters were allowed to 
"float" freely; the only exception to this was that the period was not allowed to fall below 
the minimum period already defined. Of all of the resultant converged fits, the sinusoid with 
the minimum x 2 was selected as the best fit for that pulsar. As a final check, we re-ran the 
analysis with a 4001 x 4001 grid of initial guesses; no significant changes to our fit results 
were found. 

To check that the best sinusoidal model significantly im proves the fit when compare d 
to the constant-PPA model, we applied the F-Test (see §11.4. iBevington fc Robinsonlll992l ). 
For the best-fit sinusoid for each pulsar, we calculated the F x statistic, where F x = Ax 2 /xl- 
Here, Ax 2 is the difference in x 2 between the best-fit constant PPA and the best-fit sinusoid, 
and xl is the reduced x 2 value for the best-fit sinusoid. To assess its significance, one 
calculates the probability that the sinusoidal model's improved fit would be achieved with 
random data. To accept it as formally significant, we ask that the sinusoidal fit be an 
improvement over the constant-PPA fit at the P ~ 99.7% level (approximately 3cx); i.e., we 
require that the calculated improvement in x 2 would be less than ~ 0.3% probable, given a 
random sample of data points. 

We find that nineteen of our pulsars meet or exceed this sinusoidal vs. constant-PPA F- 
Test significance level, indicating statistically significant sinusoidal PPA variations with time. 
It is important, however, to recognize that a particular pulsar's F-Test significance calculation 
is only advisory, given the uncertainty in the value of the pulse jitter variance (see §2.1.5p . 
which enters crucially as a scale factor in the reduced x 2 calculation. Nevertheless, there is 
a clear pattern of statistically significant improvement in many fits across our pulsar sample 
when a sinusoidal model is introduced. The details of the nineteen pulsars' statistically 
significant fits, along with two additional ones (see below), are listed in Table [TJ 

To highlight our assessment of the quality of the fits, we group the 81 sinusoidal fits into 
four different classes (see Tabled] for designations of all pulsars in the three highest classes). 
The top two classes have formally significant (P > 0.997) fits, while the bottom two do not. 
"Class I" pulsars are those whose sinusoidal fits we find particularly convincing. The data 
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and fits for these four pulsars are shown in Fig. [3j "Class II" consists of those pulsars where 
we find the fit less convincing, although still formally significant, for two principal reasons, 
(i) If there is a large spread in the PPA values that seems to not be accounted for by the 
sinusoid, possibly indicating non-sinusoidal variation in the PPA (see an example in Fig. HJ); 
or (ii) if the constant-PPA fit itself appears sufficient (xt ^ 1)> therefore calling into question 
the need for a sinusoid. "Class III" pulsars are those that we judge might have sinusoidal 
variation, although the fit was not a formally significant improvement over the constant-PPA 
fit. Both such pulsars' fits are displayed in Fig. "Class IV" encompasses the remaining 
60 pulsars, which have neither formal nor apparent sinusoidal tendencies. See Figure [6] for 
an example. 

3.1.1. Ultralong- Period Variations in Polarized Position Angle 

In our sinusoidal fits (cf. §3. 1 j) . none of the trial fits with long-period variations (longer 
than the total timespan of observations) is a significant improvement over the constant PPA 
fit at the 3<T level. Therefore, there is no evidence in our data of variations in PPA, including 
ramps, on these ultralong timescales. See §3.21 for a discussion of upper limits on PPA 
variations from a (constant) mean on shorter timescales. 

3.2. Upper Limits on Temporal Variations from the Mean Polarized Position 

Angle 

Table [2] presents a full summary of our PPA versus time measurements for the 81 
pulsars. The table is organized into three principal categories, (i) Pulsar Details: The first 
seven entries for a given pulsar list catalog information such as pulsar name and DM and 
RM; (ii) Data Catalog: The next four entries detail our observations, listing total timespan, 
total number of sessions, and the number rejected due to the PPA being many aA^(t) from 
the mean (see §3.2.11 for further details on the rejects); and (iii) Retained Session Results: 
These results represent important upper limits derived from the "sessions retained" after the 
editing process, and are discussed further below. 

The standard deviation of the position angle of a given pulsar, called in the "Retained 
Session Results" section of Table El is one of the most important quantities measured in 
this experiment. It represents the typical deviation of any single PPA measurement from 
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the mean position ang (For the pulsars exhibiting possible sinusoidal variations, the 
standard deviation from the sinusoid, o"^, jSine , is also presented in Tabled]). Together, these 
results represent, to our knowledge, the only published systematic measurements of time 
variations (or their upper limits) of pulsar polarized position angles. 

3.2.1. Discontinuous Changes in Polarized Position Angle 

As noted in §3.2[ there are a handful of measured Aif)(t) lying many o~Af(t) away from 
the mean value for a given pulsar. (See the "Sessions Manually Rejected" column of Table 
[2] for a listing of the pulsars having "manually" rejected sessions, and Table [3] for details on 
these sessions' properties.) There are only six such rejected sessions among all 81 pulsars, 
and never more than one per pulsar. We have examined the polarized profile in each such 
case and see no evidence for instrumental problems. These puzzling events, if real, indicate 
the passage of ionized, magnetized structures across the line of sight on timescales shorter 
than the several-month intervals between observing sessions. 

4. Discussion 

As noted above, polarized position angle variations could be caused either by phenomena 
at the pulsar (e.g., precession or magnet ospheric adjustments in the emission beam geome- 
try); or by interestellar Faraday rotation variability. The observed measurements and upper 
limits on position angle changes reported in Tables [T] - [3] can be translated (1) directly into 
changes or upper limits thereof in projected spin axis or magnetospheric orientation; or to 
(2) changes or upper limits thereof in interstellar magnetic fields via Faraday rotation. We 
now discuss each in turn. 



3 The quantity is uniquely determined from the set of measurements over all L retained epochs, 
A?p(tj); (j = 1, L). This contrasts with the situation for the individual session error bars, <7A-0(t) 01 Eq. 
[T71 which depend on the assumptions discussed in Q2. 1.51 regarding pulse jitter that are untestable with the 
current data. 
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4.1. Changes in Orientation of the Pulsar or its Magnetosphere 

4-1.1. Precession or Pulsar Magneto spheric Changes 

Much theoretical work suggests that the amplitude of the precession would be smaller 
than a degree, which would place it below the threshold of measurability for most pulsars in 
this experiment ( cf. Table [2]) , and furthermore that it would quickly dam p, at leas t in th e 
presence of rigid pinni n g or even th e slow creep of superfluid vortices [e.g., IShahaml ( 119771 ): 
Sedrakian et al.l ( 119991 ); iLinkl ( 120061 )]. Nevertheless, there are several observations of periodic 
changes in pulse shape and/or in arrival times of isolated pulsars in the literature, along with 
ours, that seem to indicate the presence of precession. In a few of the other cases, the derived 
amplitude is on the order of one to several degrees and the period is in the few-yr range, 
which are similar to our results. The best documented example is PSR B 1828-11, whose 
arrival times and pulse shape vary periodically and in a cor related fashion, indicating pre - 
cession with a period on the o r der of 1 yr and 0?3 amplitude (IStairs. Lvne. fc Shemarl 120 00). 
or ~ 3° (jLink &: Epstein 2001 ; Rezania 2003 ) or possibl y more (Akgiin. Link. &; Wasserman 



2006 



)B PSR B1642-03 also exhibits both phenomena ( Cordes fc Downs 1985 ; Blaskiewicz 



199ll ). which hav e been modeled as stemming from precession with a n amplitude of 0?5 — 0?8 



over s everal yr ([Corded Il993l ; IShabanova. Lyne. fc Uramal l200ll ). 4 ISuleymanova &; Shitov 



( 11994J ) discovered pulse shape and timing variations in PSR B2217+47, which they mod- 
eled as induced by precession with a timescale > 20 yr. PSR B1557-50 shows a several- 
year timing and dispersion measure periodicity (the latter interpreted as possible frequency- 
dependent tem poral pulse shape variations) which is m o deled as precession with an ampli- 
tude of ~ 0?01 dChukwude. Ubachukwu. fc Okekell2003l ). b'Allesandro fc McCullochl dl997h 



observed quasiperiodic timing and profile variations in PSR B0959-54 which could be caused 
by precession with a period > 2500 d and amplitude < 0?15. Quasiperiodicities on several 



(Lvne et al. 


1988; 


Scott et al. 


2003 





PPA, and timing variations in the Vela pulsar which they suggested cou ld result from free 



precesssion with a 330-d period; and iDeshpande fc Radhakrishnanl (120071 ) further elaborated 



the precession model to explain variations in X-ray morphology of the Vela Pulsar's Pulsar 



Wind Nebula. IZhu fc Xul (120061 ) have proposed that the Galactic Center Radio Transient 
GCRT J1745-3009 could be a neutron star precessing by > 15°, although an experiment like 
the current one would be insensitive to its 77-min period. Finally, the X-ray emitting iso- 
lated neutron star RX J0720. 4-3125 exhibits a wide variety of correlated periodic phenomena 



1 Lvne et al. ( 2010l ) show that pulseshape changes are correlated with spindown changes in this pulsar, 
suggesting a magnetospheric rather than precession explanation. 
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indicative of precession on a 7-8 yr timescale ( IHaberl et al.l 120061 ). 



Our analysis shows that 19 of our 81 pulsars exhibit PPA variations that are significantly 
better fit (at approximately the 3cr level) by a sinusoidal PPA function than by a constant 
PPA. Of these, we judge four pulsars to exhibit the most convincing evidence for sinusoidal 
variations in PPA. (These are our "Class I" pulsars: B0523+11, B0611+22, B0656+14, and 
B2053+21; see Fig. [3]). Our results indicate that variations, and even sinusoidal variations, 
are occasionally present in the pulsar population. As a group, these 19 pulsars display 
sinusoidal periods mostly grouped around 185 to 450 days, with a peak at 200 days and 
outliers at approximately 790, 1050, and 1250 days. The amplitude of the variation ranges 
from 1° to 8° for most of the significant fits, with a peak at 2° and an outlier at 12°. Note 
that these sinusoidal results do not include the six large PPA excursions, rejected from the 
fits in Table [2] and detailed in Table El as these large excursions do not carry the signature 
of precession even if they are real. 

In §2. 1.1b . we briefly discussed rotating magnetospheric currents as a possible source of 
periodic PPA changes. In the absence of well-developed models, we only note here that such 
phenomena might provide an alternate explanation for the observed periodic PPA variations. 



4.2. Interstellar Magnetic Field Variations 



Table [3] lists the properties of those few PPA measurements lying far from the mean. 
As noted above, these could indicate the presence of magnetized variations in the ISM, 
although they might conceivably result instead from an experimental problem. Under the 
first assumption (magnetized ISM variations causing the deviations), we calculate and list 
in Table [3] the quantities ARM(t) (Eq. EJ) and ARM(t)/RM w fAB\\(t)/(B\\) (Eq. [@. 
However, the great majority of our pulsars shows no evidence whatsoever of PPA variations 
above the noise. Under the second assumption (that the few large PPA deviations observed 
are spurious), our measured upper limits on temporal variations of position angle among 
"Retained Sessions" (the quantity of Table |2J) can be interpreted as upper limits on 
magnetized variations, again via Eqs. |3] and [7] where the quantities Aifj Faraday and ARM 
are replaced by their respective statistical analogs, and a km- Table [2] lists these derived 
upper limits on a RM and <r RM /RM (« fa B J(B\\)). 

There are occasional reports in the literature of temporal variations in interstellar mag- 
netic fields, some of which may hav e origins similar to ours. First, there are the cases of 
the Vela and Crab Nebula Pulsars ( [Hamilton et all 11985c iRankin et all 119881 ) . The varia- 
tions of RM (and also of DM) are so large and frequent that they are almost certainly 
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associated with passage of inhomogeneous nebular SNR material across the line of sight. 
This interpretation is easy to reconcile with these pulsars' relative youth and the observed 
presence of an SNR. There are other reports in the literature of RM variations toward or- 
dinary pulsars on multiy ear timescales, sometimes accompanied by DM variations as well 
(Ivan Ommen et al.lll997D. which a ppear to be caused by the passage of a cloud across the 
line of sight. IWeisberg et al.l ( 120041 ) measured many pulsar RMs and compared their results 
with earlier work wherever possible. They noted numerous cases of RMs changing signifi- 
cantly with respect to values measured ~ (5 — 20) yr earlier, while the DM tended to vary 
far less on similar timescales. 

Yet our sample of pulsars exhibits only rare evidence for RM variations over the four 
years of observation, and no evidence for sustained changes or long-term trends. We can 
suggest several possible resolutions to the disagreement between the rather common RM 
changes in the literature and our finding of essentially steady RMs in our sample: (1) 
RM variations tend to be small on the four- year timescale of our experiment, but grow on 
decade-long scales. This would have interesting implications for the fluctuation spectrum of 
interstellar plasma. However, recall that we have some sensitivity to changes on timescales 
longer than our dataspan, but no evidence was found for variations on scales beyond ~ 10 3 
d. (2) A subtle nonortho gonal emission mode comp etition process can lead to spurious 
apparent RM variations ( IRamachandran et al.l 120041 ) . which might grow with time. (3) 
One of the previously published multiepoch RM measurements was incorrect. This latter 
possibility cannot be ruled out, since RM measurements are difficult, and since the pairs of 
observations were frequently performed by different groups on different telescopes. Further 
measurements are required in order to determine the correct explanation. 



5. Conclusions 

We have studied the temporal behavior of polarized position angle in 81 pulsars over 
a ~ 4 yr period. We demonstrate that the intrinsic pulse-to-pulse jitter of polarized sub- 
pulses frequently dominates the uncertainty in position angle determination. We searched 
particularly for sinusoidal variations, discontinuous changes, and linear trends in the data. 
Nineteen pulsars have variations in PPA that were significantly better fit (at the 3a level) 
by a sinuoidal function than by a constant; such variations could be caused by precession 
or other cyclical phenomena at the pulsar. These pulsars (and especially the four pulsars 
in our "Class I") are good candidates for further study. A handful of discontinuous changes 
were observed in the sample, which could be due to the passage of the line of sight across 
an interstellar region with significantly different magnetized properties. No linear variations 
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(i.e., "ramps") of PPA across the timespan of observations were detected. For all 81 pulsars, 
we also calculate upper limits on random variations of PPA from the mean, which provide 
numerical constraints on precession and magetized interstellar variations in our sample. 

DGB, JEE, JJM, and JMW have been supported by grants from the National Sci- 
ence Foundation. JEE would like to acknowledge support from NSF grants AST 0507367, 
PHY-02 15581 (to the Center for Magnetic Self-Organization in Laboratory and Astrophys- 
ical Plasmas), PHY 0821899, and AST 0907837 to the University of Wisconsin. Arecibo 
Observatory is operated by Cornell University under cooperative agreement with the NSF 
We thank Matthew Frank and Ryan Terrien for computational assistance. 

APPENDIX 



A. Equations for Determination of a Target Pulsar's PPA Offset with Respect 

to a Calibrator Pulsar 

At an observing session at epoch t, we measure the position angle curves of target 
and calibrator pulsars as a functio n of longitude d > %l){(t> : t) and ?/> ca i(0, t), respectively. The 



high S/N position angle curves of IWeisberg et al.l (119991 ). which are coherent sums over all 



sessions, serve as templates against which to measure shifts in position angle for a given pulsar 
on a particular observing session. The position angle curve of these target and calibrator 
templates as a function of longitude 0, is ^(<f>) and ^> C ai (</>), respectively. (See Fig. |2] for 
examples of session and template position angle curves of a target and a calibrator pulsar). 

We first determine ip(t), the target pulsar's weig htecfl, longitude-averaged difference 
between PPA at an observing session at epoch t and PPA of its summed template (see Fig. 
for a target pulsar's session and template PPA curves): 



m = ^^mw-n.^ (A1) 



where the summations are over R longitude bins. 

Similarly, we also determine i[) C3 i(t), the calibrator pulsar's weighted 4 , longitude-averaged 
difference between PPA at an observing session at epoch t and PPA of its summed template 



5 The weighting factors, w((f>), account for S/N varying with longitude bin. 
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(see Fig. [2}d for a calibrator pulsar's session and template PPA curves): 

E r =l «W 

Then the quantity of interest for this experiment, Aip(t), the polarized position angle 
offset of the target pulsar with respect to a calibrator pulsar at epoch t, is: 

AV(t)=V(*)-Vw(t). (A3) 



Note then that the desired Aip(t) of a target pulsar is ultimately determined via a 
double-differencing procedure on the various PPAs - first both the target and calibrator 
pulsars' PPAs at epoch t are differenced from their templates, as in Eqs. IA1I and IA2j and 
finally the resulting target and calibrator epoch-template differences are differenced from 
one another, cf. Eq. IA31 

As pointed out in §2] (see also Eq. [T]), all of the above measured position angles contain an 
unknown instrumental origin VWtr(^)- Note that care was taken to measure both the target 
and the calibrator pulsar with an identical instrumental configuration and hence identical 
instrumental origin throughout a session at some epoch t. Therefore ipi ns tr{t) drops out of 
Eq. IA3I due to the differencing of two terms, each containing this quantity. As a result, 
temporal variations in PPA are accessible even if the absolute PPAs themselves are not. 

Further, we focus on temporal changes by subtracting a weighted mean offset Aip for 
each pulsar, before plotting the final values Aip(t), where 



A$(t) = Ai/j(t) - Ail). (A4) 

and 

£>(t)AV>(t) 



A^ = " w , (A5) 

E™w 

with the summations being over all observing sessions. The horizontal solid line placed at 
Aip = in Figs. [3] and H] then delineates the mean measured position angle in the new 
coordinate system defined by Eq. IA4I The error bar on each individual measurement is 
&Ai/>(t) > as discussed above in §2.21 and defined in Eq. ITT1 
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Table 1. Best Sinusoidal Fits to Polarized Position Angle as a Function of Time 



Pulsar 


Const. 








Sinusoidal PPA 








PSR J 


PSR B 


PPA 
























xi 


xi 


1 mm 

[d] 


p 

1 max.i 

[d] ' 


^"best.i 

[d] 


Ampl. 

a* [°1 


^ip ,sinc 

[°] 


F-Tcst 


Prob. 


Class 


J0525+1115 


B0523+11 


1.92 


0.03 


435 


6532 


1250 


12.46 


0.47 


312.51 


0.9968 




J0614+2229 


B0611+22 


1.01 


0.63 


187 


6546 


313 


2.14 


1.17 


10.73 


0.9982 




J0629+2415 


B0626+24 


0.63 


0.14 


228 


5692 


356 


2.33 


0.55 


35.25 


0.9997 


II 


J0659+1414 


B0656+14 


1.38 


0.47 


285 


7120 


1047 


2.47 


0.85 


20.24 


0.9985 




J1136+1551 


B1133+16 


1.37 


0.75 


234 


7603 


364 


2.81 


1.47 


12.97 


0.9987 




J 1607-0032 


B1604-00 


3.95 


0.48 


324 


5667 


380 


12.42 


3.14 


46.31 


0.9948 


III 


J1740+1311 


B1737+13 


1.05 


0.69 


106 


6615 


185 


2.97 


1.80 


15.94 


1.0000 




J1805+0306 


B1802+03 


1.98 


1.70 


114 


7118 


209 


2.46 


2.73 


6.95 


0.9980 




J1823+0550 


B1821+05 


1.25 


0.90 


143 


7129 
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4.62 


2.92 


10.40 


0.9995 




J1844+1454 


B1842+14 


0.60 


0.36 


158 


7118 


221 


1.67 


1.22 


14.60 


0.9999 




J1857+0057 


B1854+00 


6.08 


4.86 


150 


7129 


270 


7.67 


7.21 


7.54 


0.9974 




J1915+1009 


B1913+10 


0.41 


0.25 


190 


7123 


443 


1.69 


1.06 


12.62 


0.9993 




J1916+0951 


B1914+09 


0.23 


0.06 


265 


6634 


454 


1.92 


0.59 


27.47 


0.9993 




J1920+2650 
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2.25 


0.26 


345 


5179 


869 


15.38 


2.68 


41.51 


0.9764 


III 
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0.82 


0.44 


168 


7126 


787 


1.92 


1.23 


16.91 


0.9999 




J1954+2923 
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0.83 


0.39 


158 


7126 
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4.00 


1.84 


22.24 


1.0000 




J2022+2854 


B2020+28 


0.89 


0.47 


160 
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176 


4.21 


1.57 


18.84 


1.0000 




J2037+1942 


B2034+19 


0.75 


0.49 
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6521 
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3.69 


9.72 


0.9974 




J2055+2209 


B2053+21 


2.69 


1.42 
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3.29 
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0.9990 
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2.13 


1.24 


186 
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4.66 


11.50 


0.9980 




J2212+2933 


B2210+29 


0.70 


0.23 


201 


6045 


310 


3.65 


1.29 


26.25 


0.9998 
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Table 2. Measured Parameters and Quanities Derived Assuming Constant PPA 



Pulsar Details Data Catalog Retained Session Results 



Name DM RM 
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J0304+1932 


B0301+19 


15.74 


1 


-8.30 


0.30 


6 


3.6 


13 





13 


1.49 


0.59 


0.071 


0.6 





J0525+1115 


B0523+11 


79.34 


1 


35.00 


3.00 


7 


3.6 


6 





6 


3.73 


1.48 


0.042 


1.9 





J0528+2200 


B0525+21 
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1 


-39.60 


0.20 


8 


3.6 


16 





16 


3.74 
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1 
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1 A 

14 
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1 
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9 


1.87 
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3.46 


0.112 


3.4 


1 


J1740+1311 


B1737+13 


48.67 


1 


64.40 


1.60 


7 


3.6 


26 


1 


25 


2.44 


0.97 


0.015 


1.1 





J1805+0306 


B1802+03 


80.86 


1 


100.00 


100.00 


a 


3.9 


25 





25 


3.21 


1.27 


0.013 


2.0 


1 


J1823+0550 


B1821+05 


66.78 


1 


145.00 


10.00 


9 


3.9 


20 





20 


3.85 


1.52 


0.011 


1.2 





J1825+0004 


B1822+00 


56.62 


1 


21.00 


13.00 


7 


3.1 


5 





5 


7.40 


2.93 


0.139 


1.3 





J1841+0912 


B1839+09 


49.11 


1 


53.00 


5.00 


9 


3.9 


16 





16 


3.89 


1.54 


0.029 


2.4 


1 


J1844+1454 


B1842+14 


41.51 


1 


121.00 


8.00 


9 


3.9 


18 





18 


1.79 


0.71 


0.006 


0.6 





J1848-0123 


B1845-01 


159.53 


1 


580.00 


30.00 


10 


3.9 


10 





10 


1.28 


0.51 


0.001 


0.2 





J1851+0418 


B1848+04 


115.54 


1 


100.00 


100.00 


a 


3.6 


4 





4 


16.11 


6.38 


0.064 


3.7 





J1851+1259 


B1848+12 


70.61 


1 


158.00 


16.00 


n 


3.9 


12 





12 


5.27 


2.09 


0.013 


1.4 





J1850+1335 


B1848+13 


60.15 


1 


146.00 


8.00 


n 


3.9 


17 





17 


1.64 


0.65 


0.004 


0.5 





J1856+0113 


B1853+01 


96.74 


1 


-140.00 


30.00 


12 


3.6 


14 





14 


12.48 


4.94 


0.035 


3.4 


2 


J1857+0057 


B1854+00 


82.39 


1 


104.00 


19.00 


7 


3.9 


19 





19 


9.08 


3.59 


0.035 


6.1 


6 


J1901+0156 


B1859+01 


105.39 


1 


-122.00 


9.00 


12 


3.9 


8 





8 


8.14 


3.22 


0.026 


1.9 


1 


J1901+0331 


B1859+03 


402.08 


1 


-237.40 


1.50 


9 


3.6 


9 


1 


8 


1.05 


0.42 


0.002 


0.1 





J1903+0135 


B1900+01 


245.17 


1 


72.30 


1.00 


9 


3.6 


10 





10 


1.94 


0.77 


0.011 


0.4 





J1902+0556 


B1900+05 


177.49 


1 


-113.00 


11.00 


9 


3.6 


10 





10 


3.06 


1.21 


0.011 


0.5 





J1902+0615 


B1900+06 


502.90 


1 


100.00 


100.00 


a 


3.6 


7 





7 


9.77 


3.87 


0.039 


1.5 





J1905-0056 


B1902-01 


229.13 


1 


100.00 


100.00 


a 


3.9 


7 





7 


10.11 


4.00 


0.040 


2.2 





J1906+0641 


B1904+06 


472.80 


1 


100.00 


100.00 


a 


3.3 


9 





9 


1.85 


0.73 


0.007 


0.5 





J1909+0254 


B 1907+02 


171.73 


1 


100.00 


100.00 


a 


1.3 


5 





5 


5.79 


2.29 


0.023 


0.5 
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Table 2 — Continued 



Pulsar Details Data Catalog Retained Session Results 



Name DM RM 

















Time- 


Total 


Sessions 


Ses- 




CRM 


CRM 


Re- 


Outli 


PSR J 


PSR B 


Value 


Rcf. 


Value 


± 


Ref. 


span 


Ses- 


Manually 


sions 






\RM\ duced 


> 3 






(pc/cm 3 ) 




(rad/ 


m 2 ) 




(yr) 


sions 


Rejected Retained (deg) (rad/m 2 ) 




V 2 
A 




J1910+0358 


B1907+03 


82.93 


1 


-127.00 


7.00 


9 


3.9 


17 


o 


17 


6.53 


2.58 


0.020 


1.1 





J1909+1102 


B1907+10 


149.98 


1 


540.00 


20.00 


9 


3.9 


11 


g 


\\ 


2.13 


0.84 


0.002 


0.7 





J1912+2104 


B1910+20 


88.34 


1 


148.00 


10.00 


9 


2.5 


7 


g 




5.39 


2.13 


0.014 


0.5 





J1914+1122 


B1911+11 


100.00 


2 


100.00 


100.00 


a 


3.9 


9 


Q 


9 


7.58 


3.00 


0.030 


2.1 





J1913+1400 


B1911+13 


145.05 


1 


435.00 


30.00 


9 


3.9 


13 


n 
u 


J-O 


1 9^ 


1 29 


0.003 


0.8 





J1915+1009 


B1913+10 


241.69 


1 


431.00 


22.00 


9 


3.9 


15 


o 


15 


1.61 


0.64 


0.001 


0.4 





J1916+0951 


B1914+09 


60.95 


1 


100.00 


6.00 


9 


3.6 


10 


o 


10 


1.47 


0.58 


0.006 


0.2 





J1916+1312 


B1914+13 


237.01 


1 


280.00 


15.00 


9 


3.9 


12 


n 
u 


12 


9 Q7 
z.y i 


IIS 


0.004 


1.1 





J1917+1353 


B1915+13 


94.54 


1 


233.00 


8.00 


9 


3.9 


13 


o 


13 


2.05 


0.81 


0.003 


0.8 





J1918+1444 


B1916+14 


27.20 


1 


100.00 


100.00 


a 


3.6 


12 


I 


\\ 


4.13 


1.63 


0.016 


1.6 





J1919+0021 


B1917+00 


90.31 


1 


120.00 


7.00 


9 


2.8 


5 


o 


5 


12.19 


4.82 


0.040 


2.2 





J1920+2650 


B1918+26 


27.62 


1 


100.00 


100.00 


a 


2.8 


6 


o 


g 


7.91 


3.13 


0.031 


2.2 





J1921+1419 


B1919+14 


91.64 


1 


100.00 


100.00 


a 


3.9 


9 


n 
u 


Q 

y 


A Q'? 

^.yo 


±.yo 


0.020 


1.3 





J1921+2153 


B1919+21 


12.46 


1 


-16.50 


0.50 


9 


2.8 


7 





7 


3.53 


1.40 


0.085 


0.4 





J1922+2110 


B1920+21 


217.09 


1 


282.00 


14.00 


9 


2.5 


7 





7 


10.13 


4.01 


0.014 


2.2 





J1926+0431 


B1923+04 


102.24 


1 


0.00 


11.00 


9 


2.8 


8 





8 


2.89 


1.14 


b 


0.3 





J1926+1648 


B1924+16 


176.88 


1 


320.00 


14.00 


9 


3.6 


16 


1 


15 


1.33 


0.53 


0.002 


0.2 





J1932+2020 


B1929+20 


211.15 


1 


10.00 


6.00 


9 


2.3 


4 





4 


10.28 


4.07 


0.407 


1.4 





J1933+1304 


B1930+13 


177.90 


3 


100.00 


100.00 


a 


3.6 


5 





5 


10.09 


3.99 


0.040 


2.4 





J1932+2220 


B1930+22 


219.20 


1 


173.00 


11.00 


9 


3.9 


17 





17 


1.93 


0.76 


0.004 


0.8 





J1935+1616 


B1933+16 


158.52 


1 


-1.90 


0.40 


8 


3.6 


8 





8 


0.69 


0.27 


0.144 


0.1 





J1937+2544 


B1935+25 


53.22 


1 


100.00 


100.00 


a 


3.9 


11 





11 


1.25 


0.49 


0.005 


0.2 





J1946+1805 


B1944+17 


16.22 


1 


-28.00 


0.40 


9 


3.6 


5 





5 


3.52 


1.39 


0.050 


1.9 





J1946+2244 


B1944+22 


140.00 


4 


2.00 


20.00 


7 


2.3 


4 





4 


7.29 


2.89 


1.443 


0.4 





J1948+3540 


B1946+35 


129.07 


1 


116.00 


6.00 


9 


3.3 


9 





9 


1.51 


0.60 


0.005 


0.2 





J1952+3252 


B1951+32 


45.01 


1 


-182.00 


8.00 


7 


3.6 


11 





11 


6.55 


2.59 


0.014 


2.8 


1 


J1954+2923 


B1952+29 


7.93 


1 


-18.00 


3.00 


9 


3.9 


18 





18 


3.04 


1.20 


0.067 


0.8 





J2004+3137 


B2002+31 


234.82 


1 


30.00 


6.00 


9 


3.6 


15 





15 


2.89 


1.14 


0.038 


0.3 





J2018+2839 


B2016+28 


14.17 


1 


-34.60 


1.40 


8 


3.9 


20 





20 


3.02 


1.20 


0.035 


0.9 





J2022+2854 


B2020+28 


24.64 


1 


-74.70 


0.30 


6 


4.2 


19 





19 


2.42 


0.96 


0.013 


0.9 





J2030+2228 


B2028+22 


71.83 


1 


-192.00 


21.00 


7 


3.6 


10 





10 


4.58 


1.81 


0.009 


1.2 





J2037+1942 


B2034+19 


36.00 


5 


-97.00 


10.00 


7 


3.6 


14 





14 


5.38 


2.13 


0.022 


0.7 





J2037+3621 


B2035+36 


93.56 


1 


100.00 


100.00 


a 


3.6 


7 





7 


4.69 


1.86 


0.019 


1.4 





J2046+1540 


B2044+15 


39.84 


1 


-100.00 


5.00 


7 


3.6 


12 





12 


4.48 


1.77 


0.018 


1.0 





J2055+2209 


B2053+21 


36.36 


1 


-80.50 


3.00 


7 


3.6 


14 


1 


13 


5.45 


2.16 


0.027 


2.7 


1 


J2055+3630 


B2053+36 


97.31 


1 


-68.00 


4.00 


9 


3.6 


13 





13 


2.55 


1.01 


0.015 


1.4 


1 


J2113+2754 


B2110+27 


25.11 


1 


-37.00 


7.00 


7 


3.6 


15 





15 


4.88 


1.93 


0.052 


1.2 





J2116+1414 


B2113+14 


56.15 


1 


-25.00 


8.00 


9 


3.6 


12 





12 


7.57 


3.00 


0.120 


2.2 


1 


J2124+1407 


B2122+13 


30.12 


1 


-57.00 


8.00 


7 


3.6 


13 





13 


7.31 


2.89 


0.051 


2.1 
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Table 2 — Continued 



Pulsar Details 



Name 



DM 



RM 



Data Catalog 



Retained Session Results 













Time- 


Total 


Sessions Ses- &rm 


<?RM Re- 


Outliers 


PSR J PSR B 


Value 


Rcf. 


Value ± 


Rcf. 


span 


Ses- 


Manually sions 


\RM\ duced 


> 3ct 




(pc/cm 3 ) 




(rad/m 2 ) 




(yr) 


sions 


Rejected Retained (deg) (rad/m 2 ) 


x 2 




J2212+2933 B2210+29 


74.50 


1 


-168.00 5.00 


7 


3.3 


12 


12 2.80 1.11 


0.007 0.7 





J2305+3100 B2303+30 


49.54 


1 


-75.50 4.00 


7 


3.3 


10 


10 3.57 1.41 


0.019 1.3 





J2317+2149 B2315+21 


20.91 


1 


-37.00 3.00 


9 


3.3 


8 


8 2.79 1.10 


0.030 0.7 






a RM unknown; arbitrarily set to 100 for later calculations. 
b Undefined because RM=0. 

References. — 1: Hobbs et al. (2004b); 2: Hobbs et al. (2004a); 3: Lorimer et al. (2002); 4, Hulse & Taylor (1975); 5, Dewey 
ct al. (1988) ; 6, Manchester (1974); 7, Weisberg et al. (2004); 8, Manchester (1972); 9, Hamilton & Lyne (1987); 10, Taylor et 
al. (1993); 11, Rand & Lyne (1994); 12, Han et al. (2006). 



Table 3. Large PPA Offsets and Implied Rotation Measure Offsets 



Pulsar 



PPA Offset 



RM Offset 



PSR J 



PSR B 



(MJD) 



AV(t) 
(deg.) 



°Ay.(t) 
(deg.) 



'Ay.(t) 



ARM(t) 
(rad/m 2 ) 



AflM(t) , f AB M (t) 



<B||> 



J1543+0929 


B1541+09 


48870.4 


-31.4 


3.0 


-10.5 


-12.4 


-0.590 


J1740+1311 


B1737+13 


48678.1 


7.8 


1.9 


4.1 


3.1 


0.048 


J1901+0331 


B1859+03 


49040.5 


27.4 


3.8 


7.2 


10.9 


0.046 


J1918+1444 


B1916+14 


48675.1 


-11.0 


2.5 


-4.4 


-4.4 


-0.044 a 


J1926+1648 


B1924+16 


49159.4 


-20.8 


2.8 


-7.4 


-8.2 


-0.026 


J2055+2209 


B2053+21 


48466.8 


18.5 


3.2 


5.8 


7.3 


-0.091 



a RM unknown; arbitrarily set to 100 for this calculation. 
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L L 




(a) (b) 

Fig. 1. — Two panels from the animation starPrecession.mpg, which is available at this link . 
Each panel of the animation shows a precessing pulsar stroboscopically at times (separated 
by multiples of ~ 1 rotational period) that its beam points most directly at the observer; i.e., 
at those times when it is observable. The total angular momentum vector, L, is conserved, 
while the body symmetry axis k describes a cone about L approximately once per spin 
period. The absence of perfect synchronization between these two rates leads to a slow 
apparent precession of k at successive observer viewing times. The above two panels from 
the animation were chosen to illustrate the pulsar geometry at the extrema of its precessional 
excursion. The beam / magnetic axis is B. Note the changing orientation of the red meridian 
connecting k and B on precessional timescales. This meridian represents a fiducial magnetic 
field line projected onto the surface. In most pulsar models, linearly polarized emission will 
be oriented along (or perpendicular to) this meridian, so the precessional motion leads to a 
changing polarized position angle (PPA) on precessional timescales. 
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Fig. 2. — The polarized position angle (PPA) analysis procedure for target (a) and calibrator 
(b) pulsars. In each case, the bottom panel shows measured PPA as a function of longitude 
(f) for a single session at epoch t (thin line) and for the template (thick line), while the top 
panel shows total flux density (Stokes parameter /, top trace), linearly polarized flux density 
(L, middle trace), and circularly polarized flux density (Stokes parameter V, bottom trace) 
for the template. Note that 180° position angle jumps result entirely from the definition of 
the PPA and have no physical meaning, (a) For the target pulsar, we determine the PPA 
offset between single session and template at each longitude (|/0(</>, t) — ^>{4>)\ of Eq. IA1I) . 
and then find the weighted mean PPA offset across all longitudes, ip(t). (b) Similarly, but 
for the calibrator pulsar, we determine the PPA offset between single session and template 
at each longitude ([y> ca i(0, i) — ■0 ca i(0)] of Eq. IA2j) . and then find the weighted mean PPA 
offset across all longitudes ipcai(t)- 
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Fig. 3. — Polarized position angle (PPA) as a function of time for the four "Class I" (best) 
fits. For each pulsar, the left panel presents the data and best sinusoidal fit as a function of 
epoch, while the right panel displays the data and fit after folding modulo the best period. 
(The probability that the sinusoid gives a significant improvement over a constant PPA is 
listed in the top right of this panel.) The error bars represent the quantity ct^m, calculated 
as in Eq. [TTl Normal sessions are depicted with a filled diamond, while those which required 
an artificial composite calibrator are marked with an open square. 




Figure [3] (continued) 



- 29 - 



3= 

o 

03 

C 
< 



M 
O 
Q_ 



60 



40 



i 20 
o 

CL 



-20 





1 1 




1 .... 1 

PSR B1 854+00^ 










< 


















1 






i 



500 



1000 1500 
Days from Jan 1988 



2000 



3= 

o 

< 



o 

Q_ 

c 
o 



o 

Q. 




0.0 



0.2 



0.4 0.6 
Phase 



0.8 



1.0 



Fig. 4. — As in Figure El but an example of a sinusoidal fit that we judged to be "Class II." 
The fit had similar statistical significance as Class I fits, but we nevertheless judged it to be 
less convincing. Despite the formal success of the fit, the PPA of numerous points lies outside 
the amplitude range of the fitted sinusoid and another point appears to be entirely out of 
phase with it. See Fig. |3] caption for further explanations of the plot labels and symbols. 
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Fig. 5. — As in Figure |3j but for the two pulsars in "Class III." Pulsar B1604-00 has a 
significance-of-fit that is near our cutoff (the probability is 0.995). Given the uncertainties 
in this experiment, and the degree of agreement between the fit and observations, we include 
it as another possible example of precession. This class also include B1918+26; the prefered 
sinusoid fit here is of lower significance (owing partly to one less data point than B1604-00), 
but the correspondence between the fit and data points again lead us to flag it as perhaps 
showing precession. 
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Fig. 6. — As in Figure [31 but for a pulsar in "Class IV": in this case, the sinusoidal fit is not 
a significant improvement over a constant PPA. 



